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The neutralino-proton cross section is examined for supergravity models with 
R-parity invariance with universal and non-universal soft breaking. The re- 
gion of parameter space that dark matter detectors are currently (or will be 
shortly) sensitive i.e. (0.1 — 10) x 10~ 6 pb, is examined. For universal soft 
breaking (mSUGRA), detectors with sensitivity o~^o_ p > 1 x 10~ 6 pb will be 

able to sample parts of the parameter space for tan (3 ~ 25. Current relic den- 
sity bounds restrict m,-o < 120 GeV for the maximum cross sections, which 
is below where astronomical uncertainties about the Milky Way are relevant. 
Nonuniversal soft breaking models can allow much larger cross sections and 



can sample the parameter space for tan (3 ~ 4. In such models, mo can be 
quite large reducing the tension between proton decay bounds and dark mat- 
On. ter analysis. We note the existance of two new domains where coannihilation 
r~| ! effects can enter, i.e. for mSUGRA at large tan /3, and for nonuniversal models 
with small tan (5. 

X 

c3 ■ I. INTRODUCTION 

Supersymmetric models with R-parity invariance generally predict the existance of dark 
matter relics from the Big Bang. Experimental bounds on exotic isotopes strongly imply that 
the lightest supersymmetric particle (LSP), which is absolutely stable by R-parity invariance, 
must be electrically neutral and weakly interacting. The minimal supersymmetric model 
(MSSM) then has two possible candidates, the lightest neutralino (x?) and the sneutrino 
(z>). In gravity mediated supergravity (SUGRA) grand unified models (GUTs) [1], the 
allowed region in the supersymmetry (SUSY) parameter space where the v is the LSP is 
generally small. The absence of the decay Z —>■ v + v at LEP implies > 45 GeV and LEP 
bounds on the light Higgs is sufficient to eliminate the v as the LSP for the minimal SUGRA 
model (mSUGRA) [2] . If one further includes cosmological constraints, the sneutrino is also 
excluded for the general MSSM [3]. Thus for these models, the \\ is the unique candidate 
for cold dark matter (CDM). It is an appealing feature then of SUGRA models that the 
predicted amount of relic density of neutralino CDM is consistent with what is observed 
astronomically for a significant part of the SUSY parameter space. 
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Since the initial observation that the Xi represented a possible CDM particle [4] and 
the subsequent suggestion that local Xi hi the Milky Way might be observed by terrestrial 
detectors [5], there has been a great deal of theoretical analysis and experimental activity 
concerning the detection of local CDM particles. Recent theoretical calculations in Refs. 
[6-25] have made use of a number of different SUSY models. Thus Refs. [6-10] assume the 
MSSM model, and calculations in Refs. [11-19] are performed using mSUGRA GUT models 
(with universal soft breaking at the GUT scale M G = 2 x 10 16 GeV). Refs. [20,21,24] allow for 
nonuniversal soft breaking in the Higgs sector and Refs. [22,23,25] include also nonuniversal 
effects in the third generation. In addition, different authors limit the parameter space 
differently. 

The neutralino-nucleus scattering amplitude contains spin independent and spin depen- 
dent parts. However for detectors with heavy nuclei, the spin independent part dominates. 
In these, the neutron and proton scattering amplitudes are approximately equal, which al- 
lows one to extract from the data the spin independent neutralino-proton cross section cr^o_ p . 
The sensitivity of current experiments (e.g., DAMA, CDMS) is approximately (1 — 10) x 10~ 6 
pb for cr^o_ p [26], and perhaps a factor of 10 improvement may be expected in the near fu- 
ture. It is the purpose of this paper to examine what part of the SUSY parameter space can 
be tested with such a sensitivity. We do this by examining the maximum theoretical cross 
section that can lie in the domain 

0.1 x 10~ 6 pb < cr^o _ p < 10 x 10" 6 pb (1) 

as one varies SUSY parameters (e.g., tan/5, m^o). Our calculations are done within the 
framework of SUGRA GUT models with non-universal soft breaking allowed in both the 
Higgs and third generation squark and slepton sectors. (As discussed in Ref. [22] and will 
be seen below, it is necessary to include both Higgs and third generation nonuniversalities 
as the two can have constructive interference.) We also update earlier analyses by including 
the latest LEP bounds on the light chargino (xf) an d light Higgs (h) mass (m^± > 94 GeV, 
m h > 95 GeV) and include the b — > s + 7 and Tevatron constraints. 

In calculating cr^o_ p , we restrict the SUSY parameter space to be consistent with the 
astronomical estimates of the amount of relic CDM. This is conventionally measured by the 
quantity ^cdm = Pcdm/Pc where pcdm is the mean CDM mass density, and p c = 3Hq/8hGn 
(H = Hubble constant parameterized by H = (100 km s -1 Mpc" 1 )^, Gn = Newton 
constant). Recent measurements of Q m = p m j ' p c (p m is the matter density), H , the Cosmic 
Microwave Background (CMB), supernovae data, etc., indicate that ^cdm is smaller then 
previously thought. We will see below that both the upper and lower bounds on ^i^oh 2 
strongly affect the predicted values of cr^o_ p . 

In our analysis below we use the one loop renormalization group equations (RGE) [27] 
from M G to the t-quark mass m t = 175 GeV and impose the radiative breaking constraint 
at the electroweak scale. We start with a set of parameters at M G , integrating out the heavy 
particles at each threshold, and iterate until a consistent spectrum is obtained. One loop 
corrections are included in diagonalizing the Higgs mass matrix, and L-R mixing is included 
in the sfermion mass matrices so that large tan f3 may be treated. Naturalness constraints, 
that the gluino mass obey m- g <\ TeV, the scalar mass m < 1 TeV and |A /m | < 5 are 
imposed. Gaugino masses are assumed universal at Mq, and possible CP violating phases 
are set to zero. Thus we do not treat here D-brane models [28,29] (which will be discussed 
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in a subsequent paper). The SUSY mass spectrum is also constrained so that coannihilation 
effects are negligible. (We find, in fact, that this is a significant constraint with nonuniversal 
soft breaking even for low tan/3.) We examine tan/5 in the range 2 < tan/5 < 50, and 
include leading order (LO) corrections to the 5 ^ s + 7 decay and correct approximately 
for NLO effects [30,17]. We require that the theoretical branching ratio lie in the range 
1.9 x 10~ 4 < BK(B — > Xs'-f) < 4.5 x 10~ 4 , and use one loop corrections to the 6-quark mass 
so that rrib takes on its experimental value mb(mb) = (4.1 — 4.5) GeV [31]. (See Appendix. 
The loop correction is significant for large tan /3 and stems from the part of the Lagrangian 
given by -^\ b b L b ] R H^ + h.c.) We do not assume any particular GUT group constraints 
and do not impose b — r Yukawa unification (since the latter is sensitive to possible unknown 
GUT physics). 

In Sec. 2 we discuss the range of the astrophysical parameters that enter into the relic 
density analysis, and also the uncertainties of the quark content of the proton which affect 
our calculation of cr^o_ p . In Sec. 3 we examine the mSUGRA model where it is seen that 
cfyO-p > 1 x 10~ 6 pb (the current experimental sensitivity) requires tan/5 to be quite large, 
though this is somewhat relaxed in the domain 0.1 x 10 -6 pb < cr^o_ p < 1 x 10~ 6 pb. In 
Sec. 4 we discuss the nonuniversal models, and see that here Eq. (1) can be satisfied for 
relatively small tan/3 and large tuq, and also that cr^o_ p sustains for large m^o. The SUSY 
mass spectrum expected for our domain of cross sections is also examined. Conclusions 
are summarized in Sec. 5. A brief qualitative discussion is also given there of the effect of 
these results on proton decay since the above nonuniversal results appear to releave some of 
the tension previously noted [23] between c^-p in the range of Eq. (1) and current Super 
Kamiokande proton lifetime bounds [32]. 



II. ASTRONOMICAL AND QUARK PARAMETERS 

The basic experimental quantity that controls the SUSY analysis of relic density is Q^oh 2 . 
Recent measurements at the Hubble Space Telescope using a number of different techniques 
has led to a combined average of [33] 

H = (71 ± 3 ± 7) km s" 1 Mpc" 1 . (2) 

There is now sufficient data on the CMB anisotropies to show that VL tot ~ 1 and fl m small is 
strongly favored [34]. Measurements on clusters of galaxies yield Q m < 0.32 ± 0.05 [35,34], 
and these results are consistent with the supernovae data [36] . An analysis of combined data 
(excluding microlensing) yields [37] VL m = 0.23 ±0.08. In view of possible systematic errors, 
we assume here VL m = 0.3 ± 0.1, and since the baryonic content is Vl B = 0.05, we take 

Q^o = 0.25 ±0.10 . (3) 

Combining errors in quadrature then yields Vt^oti 2 = 0.126 ± 0.052. In the following we will 
restrict the range of by what is approximately 2 std. around the mean: 

0.02 < Vt^h 2 < 0.25 . (4) 
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(As pointed out in Ref. [17], the lower bound is the minimum amount of DM to account 
for the rotation curves of spiral galaxies.) Future measurements by the MAP and Planck 
sattelites will greatly reduce these errors. 

The fundamental SUSY Lagrangian allows one to calculate the neutralino-quark scatter- 
ing amplitude. To obtain the Xi ~ P cross section one needs to know in addition, the quark 
content of the proton. In the notation of Ref. [10], the two parameters that enter sensitively 
are 

;=^, (5) 

m p 



and 



j = (p\m s ss\p) 



m p 



where a^N is the (ir — N) cr-term and is given by 



/ can be written as [10] 



where 



and 



v-kN = ^(rn u + m d )(p\uu + dd\p) . (7) 



2 m p 



r — T7 r , (9) 

^(m u + m d ) 



_ (P\SS\P) _ 1 (Tp ^ 

\{p\uu + dd\p) a nN 

The quark mass ratios are fairly well known, and we use in the following r = 24.4 ± 1.5 [38]. 
Recently, the uncertainties in a n N and o"o have been analyzed in Ref. [10]. They find 

40 MeV ~ ct^n ~ 65 MeV, 30 MeV ~ a ~ 40 MeV . (11) 
In the following we will consider two possible choices for o- wN and a : 

Set 1: a nN = 40 MeV, a = 30 MeV; / = 0.0480, / = 0.195 . , , 

Set 2: a nN = 65 MeV, a = 30 MeV; / = 0.0693, / = 0.455 . ^ ' 

Set 1 corresponds approximately to the original analysis of Ref. [39] (and is the most con- 
servative possibility) while Set 2 is similar to the Set 2 of Ref. [10]. (Set 3 of Ref. [10] gives 
considerably larger cross sections.) In the following, we will use mostly Set 2 in showing our 
result, but we will exhibit the difference between Set 1 and Set 2 in one case to illustrate 
some of the uncertainties that exist. 
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III. THE MSUGRA MODEL 



We begin our analysis by examining the minimal SUGRA model which depends on 
four parameters and the sign of the Higgs mixing parameter \i. A convenient choice of 
parameters is m (the universal scalar mass at Mq), mi/ 2 (the universal gaugino mass at 
Mq), A q (the cubic soft breaking mass at Mq) and tan/5 = (H 2 )/(H 1 ) (where (#1,2) gives 
rise to (down, up) quark masses). It is convenient sometimes to replace rrii/2 by the gluino 
mass nig = (a3/«c) m i/2 (®g — 1/24 is the GUT scale gauge coupling constant) or m^o 
which also scales with mi/ 2 . Our sign convention on the /i parameter is defined by the 
quadratic term in the superpotential 

= iiH x H 2 = n(H°H° - H^H+). (13) 

(With this convention, the b — > s + 7 constraint eliminates mostly /i > 0.) 

In calculating cr^o_ p , one must impose the relic density constraint of Eq. (4). This is 
governed by the Boltzmann equation describing Xi annihilation in the early universe [40]: 

dn^o R 

~dt~ + 3 R Ux " = ( tT ^ v ^i)( n x° 1 ~ ( 14 ) 

where n^o is the number density of Xi, n eq its equilibrium value, a ann is the annihilation cross 
section, v re \ the relative velocity, and ( ) means thermal average. The diagrams governing 
a ann are shown in Fig. 1. The final relic density is given by 

ivf 
K'dx( 




l^k' = 2.48 x 10-" ^ ,,,J_ ... ; (15) 



where Xf = kTf/m^o ~ 1/20, T f is the freezeout temperature, Nf is the number of degrees 
of freedom at freezeout, and (T^o/T 7 ) 3 is the reheating factor. 

The relic density decreases with increasing annihilation cross section, and in order to 
understand some of the results obtained below, we first discuss which parameters control 
a ann . From Fig. 1 one expects a ann to fall with increasing m^o and also increasing m (since 
m 2 j increases with ml). However, if 2m^o is near m h , m H or (but lies below), the s- 
channel pole gives rise to a large amount of annihilation (which can reduce Q^oti 2 below the 
allowed minimum) and due to the thermal averaging, this effect can be significant when 2m^o 
is less than the Higgs mass and within five times the Higgs width of the Higgs mass [41,42]. 
The LEP data, has eliminated most of this effect for the light Higgs. However, we will see 
that since H and A become light at large tan/3, effects of this type become significant in 
that regime. Further, if one of the sleptons or squarks becomes light i.e. ~ 100 GeV, the 
t-channel annihilation will drive Q^oti 2 down [43]. This effect can again become significant 
at large tan f3 where large L-R mixing in the sfermion mass matrices reduces r?r~. 

We turn next to &^o_ p , which is governed by the diagrams of Fig. 2. We see here that the 
cross section can become large for light (first generation squarks) and light Higgs bosons. 
These regions of parameter space are just the ones that reduce VL^oh 2 , and so there can be 
a bound produced on cr^o_ p so that Q^oti 2 does not fall below its minimum. 

In order now to see the sensitivity of current detectors to mSUGRA we plot in Fig. 3 (for 
Set 2 parameters of Eq. (12)), the maximum value of cr^o_ p for tan (5 =20, 30, 40 and 50 as a 



5 



function of m^o (obtained by allowing all other parameters to vary subject to the constraints 



experiment is thus sensitive to mSUGRA for tan/? ~ 25 (i.e. c^o-p ~ 1-0 x 10~ 6 pb). We 
note that the fall off is less severe for tan (3 = 50, since at this high value of tan (3 the H and 
A Higgs become relatively light enhancing the Xi~ V cross section. This can be seen in Fig. 4 
where we have plotted m H for tan/3 = 30 and tan/3 = 50. We note also the importance of 
including the loop corrections to m b for large tan (3 (e.g. tan (3 = 50) to obtain the correct 
results here. 

Fig. 5 shows the sensitivity of the calculations to the choice of particle physics parameters. 
Set 1 gives cross sections about a factor of 2 smaller than Set 2. In the following, we will 
use Set 2 in all our analysis. 

Fig. 6 shows Q^oti 2 vs m^o for tan/3 = 30. We see, as expected, Q^oti 2 is an increasing 
function of m^o (since o ann is a decreasing function). The upper bound on Q^oh 2 then implies 

an upper bound of the neutralino mass of m^o = 120 GeV (i.e. m~ g ~ 900 GeV) as has been 
discussed previously [12-14]. Note that one can obtain cross sections within the DAMA 
sensitivity range without going to the edges of the parameter space in fi^oh 2 . Also, these 
cross sections all fall below the current experimental sensitivity before the uncertainties in 
the Milky Way astronomical parameters discussed in [8,9] become important. 

Figs. (7-9) exhibit the particle spectrum expected for the example of tan/3 = 30, when 
c^o_ p takes on its maximum value. Fig. 7 shows that the ci-squark is quite heavy (m is 
large). This arises from our constraint mf R — m^o > 25 GeV to prevent coannihilation 
effects from occurring. Thus the large tan/3 being considered here reduces (due to 
L-R mixing) and uiq must be increased to prevent it from becoming degenerate with the xi- 
These coannihilation effects are thus different from the ones that can occur at low tan (3 [20], 
since they occur at low m^o where cr^o-p is large enough to fall within the range of Eq. (1). 
(They will be discussed elsewhere.) Fig. 8 shows the light Higgs mass, which is relatively 
heavy due to the fact that m is large. Fig. 9 shows m-± vs m^o for tan/3 = 30. One sees 
that scaling is obeyed [44], i.e. m~± = 2m^o since ji is relatively large (/i 2 /M| 3> 1). 



Nonuniversal soft breaking can arise in SUGRA models if in the Kahler potential, the 
interactions between the fields of the hidden sector (that break supersymmetry) and the 
physical sector are not universal. Nonuniversalities allow for a remarkable increase in the 
neutralino-proton cross section. 

In order to suppress flavor changing neutral currents (FCNC), we will assume the first 
two generations of squarks and sleptons are universal at Mq (with soft breaking mass m ) 
but allow for nonuniversalities in the Higgs and third generation. Thus we parameterize the 
soft breaking mass at Mq as follows: 




IV. NONUNIVERSAL MODELS 




(16) 
(17) 
(18) 




m 2 (l + 5 3 ); m 2 R = 
m l = 



mg(l + 5 4 ); m 2 R = ml(l + <5 5 ); 
ml(l + 5 6 ); 171^ = 1711(1 + 57); 



6 



where m is the universal mass of the first two generations, qi = (iz,,&i,); II = 
ur = t R ; e R = tr etc. We take here the bounds 

-1 < St < 1 (19) 

An alternate way of satisfying the FCNC constraint is to make the first two generations 
very heavy, and only the third generation light. This is essentially included in the above 
parameterization by making m large, and taking the Si sufficiently close to -1. 

One of the important parameters effected by the nonuniversal soft breaking masses is /j, 2 , 
which is determined by the radiative breaking condition. While the RGE must be solved 
numerically, an analytic expression can be obtained for low and intermediate tan (3 [22] : 



t 2 

.2 1 



n = 



f l-3£>o , M , \ l ~ D » ,x i x \ + 6l ~ 



t 2 - 

+universal parts + loop corrections (20) 



where t = tan (3 and 



c - 1 - <2£b> 2 < 21 > 



A similar expression holds for large tan/3 in the 5*0(10) limit so Eq. (20) gives a qualita- 
tive picture of the effects of nonuniversalities in general (a result borne out from detailed 
numerical calculations). 

We see first that in general D Q is small, i.e. for m t = 175 GeV, D < 0.2, and hence the 
squark nonuniversality, S3 and S4, produce comparable size effects as the Higgs nonuniver- 
salities Si and S 2 , so that both must be included for a full treatment [22]. Second, one can 
choose the signs of Si such that either /i 2 is reduced or /i 2 is increased. The significance of 
this is that in general, the Xi is a mixture of Higgsino and gaugino pieces 

Xi = aW 3 + pB + 7^ + SH 2 (22) 

Now the spin independent part of Xi ~ 1 scattering depends on interference between the 
gaugino and Higgsino parts of x± [39] (it would vanish for pure gaugino or pure Higgsino) and 
this interference increases if \j? is decreased (increasing cr^o_ q ) and decreases if jj 2 is increased 
(decreasing cr^o_ q ) [15]. Thus there are regions in the parameter space of nonuniversal models 
where cr^o_ p is significantly increased compared to the universal case. 

The above effect can be seen in Fig. 10 where the maximum cr^o-p are plotted for tan (3 = 
7 for the nonuniversal and universal cases. We see that nonuniversalities can increase cr^o_„ 
by a factor of ~ 10. Fig. 11 plots the nonuniversal curves for tan/3 = 3,5, and 7. One sees 
here that with nonuniversal soft breaking, the current DAMA sensitivity requires tan j3 ~ 4 
(compared to tan/3 ~ 25 in the universal case). For larger tan/3 one can get very large 
nonuniversal cross sections. Fig. 12 shows the maximum cr^o_ p for tan/3 = 15, which already 
lies in the region excluded by CDMS and DAMA. 

For GUT groups containing an SU(5) subgroup (such as SU(5), S'O(IO), SU(6) etc.) 
with matter in the usual 10 + 5 representations, the Si of Eqs. (17,18) obey 

S 3 = S 4 = 5 5 = 5 W ; 5 e = 5 7 = <5§ (23) 
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We consider this case in more detail (where it is assumed that the gauge group breaks to 
the Standard Model at Mq)- Fig. 13 shows Q^oh 2 when cr^o_ p takes on its maximum value 
for the characteristic example of tan/3 = 7. One sees that Q^oh 2 is generally small since one 
has 5 W < to obtain the maximum cr^o_ p (reducing /i 2 in Eq. (20) and hence increasing the 
cross section). This however reduces m fR (from Eq. (17)) increasing the annihilation rate as 
in the discussion of Fig. 1. (If the S77(5)-type constraint were relaxed and 5 5 left arbitrary, 
VL^nti 1 could be increased. For example, <5 5 = produces m 50% increase in Vt^oh 2 .) The 

further fall off of ^j^h 2 for m^o ~ 110 GeV arises from the fact that m H ~ 300 GeV, and the 
nearness of the s-channel pole of Fig. 1 increases the early universe annihilation. This 
can be seen explicitly in Fig. 14 where 2m^o is close to m# when m^o ~ 110 GeV. Fig. 15 
shows that the light Higgs for this case is quite light lying just above the LEP2 bounds. 
Particularly interesting is that the first two generations of squarks, however, are relatively 
heavy. This is shown in Fig. 16 for the <i-squark. The reason for this can be seen from Eq. 
(20) where since 5$ = 64 = 5 10 < (to lower yU 2 and hence increase <y^o_, p ) the nonuniversal 
terms produce a net negative m 2 , contribution to /i 2 , the lowering of ji 2 being enhanced, 
then, the larger m is. Thus it is possible to get heavy squarks in the first two generations 
at low tan /3, which may have implications with respect to proton decay as discussed in the 
next section. 



V. CONCLUSIONS 

If the dark matter of the Milky Way is indeed mainly neutralinos, then current detectors 
are now sensitive to interesting parts of the SUSY parameter space. Thus either discovery 
(or lack of discovery) will determine (or eliminate) parts of the parameter space, and this 
analysis is complementary to what one may learn from accelerator experiments. 

To examine what parts of the parameter can be tested with current detectors or in the 
near future, we have considered cr^o_ p , the Xi~~P cross section, in the range 0.1 x 10~ 6 pb < 
a xl'P < 10 x 10~ 6 pb, and have plotted the maximum theoretical cross section for different 
SUGRA models. There is a major difference between the universal and nonuniversal soft 
breaking models. Thus the current DAMA experiment (with sensitivity of o"y0_ p ~ 1 x 10- b 

pb) is sensitive to tan/3 ~ 25 for universal soft breaking (Fig. 3) while it is sensitive to 
tan/3 ~ 4 for the nonuniversal model (Fig. 11). Thus while dark matter cross sections 
increase with tan/3 and hence detectors are more sensitive at higher tan/3, it is possible 
for current detectors to probe part of the low tan/3 parameter space for the nonuniversal 
models. 

For the mSUGRA model, we find that Q^oh 2 monotonically increases with m^o from 
the minimum to the maximum bounds of Eq. (4) (Fig. 6), leading to the upper bound 
m^o < 120 GeV (nig ~ 900 GeV) which is below where astronomical uncertainties about 
the Milky Way [8,9] become significant. In general /z 2 is large (i.e. /i 2 /M§ ^> 1) leading to 
the usual gaugino scaling relations e.g. Fig. 9, and the Higgs mass is relatively heavy (Fig. 8). 
At the very largest tan/3, e.g. tan/3 = 50, the loop corrections to A5 at the electroweak scale 
become very large (see Appendix), requiring that the b- Yukawa coupling, be adjusted so 
that one obtains the experimental 6-quark mass [31]. 
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For the nonuniversal model, significantly increased Xi~P cross sections can be obtained 
by choosing 5 3j 4 < and 62 > in Eq. (20). This reduces /i 2 , increasing the Higgsino 
content of the Xij an d hence increasing the Higgsino-gaugino interference which enters in 
cr^o-p- (In the S77(5)-like models, this generally leads to a light f R and hence a relatively 
low (Fig. 13)). In this case the maximum cross sections arise with /i 2 relatively small, 

and so scaling no longer holds accurately, and the light Higgs lies close to the LEP2 bounds 
(Fig. 15). 

While coannihilation effects have not been treated in this analysis, we have noted two 
regions where such effects can occur. In mSUGRA models, due to the fact that tan /3 must 
be large to obtain u^o_ p in the range of Eq. (1), L-R mixing reduces m,f R making the tr 
near degenerate with the Xv I n t ne nonuniversal case, where tan/3 is small or moderate, 
large CyV-p are obtained by lowering /i 2 which makes the xt nearly degenerate with the Xi- 
Both these domains of coannihilation are different that previously treated [19], and they 
inhabit regions of parameter space with cr^o_ p within the reach of current detectors. We 
have prevented coannihilation here from becoming significant by imposing the constraints 
m fR — m^o, m~± — m^o > 25 GeV. Further study is required to see what occurs when these 
constraints are removed. 

It has for sometime been realized that tension exist in GUT theories that simultaneously 
allow for dark matter and proton decay [23,45]. Thus for SU (5)-type models, minimal 
SUGRA GUT proton decay proceeds through the H 3 , the superheavy Higgsino color triplet 
components of the Higgs 5 and 5 representations. The basic diagram is shown in Fig. 17, 
showing that the decay rate scales approximately by 

^^'^(fe)' (24) 

where M 3 = 0(M G ) is the H 3 mass. In mSUGRA models, scaling is generally a good 
approximation and m^± = 2m^o. Hence proton stability requires small m^o, large trig and 
small tan/5. We have seen, however, that if dark matter exists with the sensitivity of the 
current DAMA experiment, while moderately heavy squark masses could exist in mSUGRA 
(Fig. 7), tan/3 would have to be quite large i.e. tan/3 ~ 25, which would be sufficient to 
violate the current Super Kamiokande bounds on the proton lifetime [32]. However, this 
tension is releaved for the nonuniversal SUGRA GUT models. Thus we saw in these cases, 
one could have a cr^o_ p in the range of the DAMA experiment for small tan /3, i.e. tan /3 ~ 4, 
and further such large cross sections also implied large squark masses, Fig. 16. This would 
be expected to remove any disagreement between a large cr^o_ p and a small proton decay 
rate. 

Finally we mention that in this paper we have plotted the maximum xi — P cross sections 
for each tan/3 and m^o. Of course nature may not chose SUSY parameters such that cr^o_ p 
takes on its maximum value. However, by looking at the maximum cr^o_ p we are able to 
see in a given model whether detection of dark matter at current detector sensitivities is 
consistent with the predictions of the theoretical model. 
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VII. APPENDIX 

The 6-quark coupling to the down type Higgs field which gives rise to tree level bottom 
mass is described by 

L bbH = \ b b L b R H° + h.c. (25) 

There also exists a term in the Lagrangian where the bottom squarks are coupled to the up 
type neutral Higgs an d is given by: 

L~ bbH = -X bf i*b L b R Hf + h.c. (26) 

The above interaction can give rise to a one loop contribution to the tree level bottom 
mass [46]. We do the analysis in the mass insertion approximation which produces errors 
of less than 10% in m b for the relevant parts of the parameter space. The loop diagram 
arising from the above interaction, shown in Fig. 18a, involves gluino, squark fields, a s and 
tan (3 and hence can be large for large tan f3. There also exists another one loop contribution 
which involves the stop quarks and the chargino. This loop, shown in Fig. 18b, depends 
on and contributes less than the gluino loop. The net 6-quark mass generated from the 
above contributions is m b + 5m b , where 

6m b = \ b v 1 KtfwP; v 1 = (H®) (27) 



where 

_ = "Htl^Na^NS] (29) 

(a — b){b — c)(a — c) 

rrig is the gluino mass, m bh are the left and right handed sbottom masses and mi L R are 
the left and right handed stop masses. 

The correction K is evaluated at the electroweak scale which we take here to m t (the 
endpoint of running the RGE down from Mq). Using the RGE for X b , we then determine 
X b (m t ) so that the total 6-quark mass, m b = X b vi + 5m b , agrees with the experimental value 
of m b {m b ) [32] at the b scale. This produces a significant change in X b for large tan/3. 
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FIGURES 
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FIG. 1. Diagrams for early universe annihilation of Xi through Higgs (h, H, A) and Z poles 
and squark and slepton (/) poles. 
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FIG. 2. Diagrams for neutralino quark scattering. 
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FIG. 3. The maximum value of <J^o_ ri vs. m^o for tan/3 = 20, 30, 40, and 50 for Set 2 
parameters of Eq. (12). 
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FIG. 4. rriH vs. m^o when cr^o_ p takes on its maximum value. Top curve is for tan j3 = 30, 
bottom curve for tan/3 = 50. 
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FIG. 5. Maximum cr^o_„ vs. 
parameters (dashed). See Eq. (12). 
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FIG. 7. m t vs. m^o for tan (3 = 30 when c^o takes on its maximum value. 
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FIG. 9. m~± vs. TTL-o for tan/3 = 30 when <r,~o_„ takes on its maximum value. 
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FIG. 10. Maximum value of cr^o_ p vs. m^o for tan/3 = 7 for nonuniversal soft breaking (upper 
curve) and universal soft breaking (lower curve). 
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FIG. 12. Maximum value of cr^o_ p vs. m^o for tan/3 = 15 for nonuniversal soft breaking. 
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FIG. 13. Q-^oh 2 vs. m^o for maximum cr^o_ p for tan/3 = 7, nonuniversal soft breaking. 
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FIG. 16. rag vs. m^o for maximum P%o_ p for tan (3 = 7, nonuniversal soft breaking. 
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FIG. 17. Example of p-decay diagram. The major contribution comes from the second gen- 
eration loop (i=2), the third generation contributing ~ 30% correction with arbitrary relative 
phase. 
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FIG. 18. One loop correction to 6-mass. 
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